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Abstract 

\^ • The arrival time coherence of particles in the Ultra High Energy Air Showers where the center 

I of mass energy of the interaction is of the order of lO^^e^, puts strict constraint on the propagation 

■ of particles in a hypothetical extra-dimension. We first argue that at such high energies bulk modes 

and massive KK-modes can be produced abundantly and in many models their phase space volume 
^ . is larger than confined modes. Then, we study the minimum propagation time in one and two-brane 

models and show that a large part of the parameter space of these models arc ruled out unless the 
confinement of fields is proteced by symmetries up to energies not accessible even to the high energy 
\^ ' tail of Ultar High Energy Cosmic Rays (UHECRs). As a by-product we confirm the result obtained 

in some previous works about the close relation between a small Cosmological Constant and the 
hierarchy problem. 

m ■ 
> 

^ ■ 1 Introduction 
O 

. for unification of Gravity with Electromagnetism [1], space-times with more than 4 dimensions have 
I been the hope of physicists to solve problems of High Energy Particle Physics. Last ideas in these 
Q^' series are suggestions by N. Arkani-Hamed, et al. [2] and by L. Randall and R. Sundrum [3] for 



Oh! 



Ever since the proposition by Th. Kaluza and O. Klein in 1920s to use a 5-dimensional space-time 



using large extra dimensions and localized matter to solve mass hierarchy problem inspired by some 



D , previous works of V. Rubakov and M. Shaposhinkov [4] on domain walls in higher dimensional spaces 
' and P. Horava and E. Witten [5] on M-theory models with Compactification in spaces with D-brane 
boundaries. 

^ \ In the first proposals only gravity could propagate in the bulk. It has been however found that the 
' total localization of all fields except graviton on 3-branes is not realistic. In fact brane solution are 
cosmologically unstable and at least one scalar bulk field (radion) [6] [7] is necessary to stabilize the 
distance between branes. In some brane models inflaton [8] also has to propagate to the bulk to make 
inflation with necessary properties. A deeper insight to the propagation of gravitational waves and 
massive particles with bulk modes in models with infinite bulk has illustrated that even the warping 
of the bulk can not stop their escape from branes [10] [11]. 

Most of localization mechanisms are evolutionary i.e. based on special configuration of matter fields 
with localized properties like topological defect solutions which can arise during phase transitions in 
the Early Universe [12]. In these models the real dimensionality of the space-time is larger than 4 but 
our Universe is confined to a domain wall (3-brane) where fields specially at low energies (with respect 
to quantum gravity scale) are geometrically or gravitationally localized. However, geometrical settings 
like a warped metric are not always enough to confine fields on the branes. It has been shown that in 
spaces with > 3 extra-dimensions gravity can not be geometrically confined to a 3-brane and p-form 
fields in the bulk must be added to stabilize the brane (defect) [13]. In 5-dim. models gravity and scalar 
fields can be localized on the brane with negative tension [14] (or the brane with smallest value of warp 
factor in 2-brane models of [15] [16]). Warp geometry can localize fermions only on the positive tension 
brane (which can not solve the hierarchy problem). Vector fields can not be geometrically confined. 
Localization of gauge fields and fermions on the negative tension brane is achievable through special 
particle physics setups [14] [17]. The localization scale however is considered to be not much higher 



than warping scale, otherwise a new hierarchy can appear [18]. At higher energies one expects that 
symmetry restoration (e.g. chiral symmetry of fermions) leads to escape of particles from the brane. 

Even when a warped geometry is enough to confine fields on the brane, the wave function of the zero 
mode can penetrate to the bulk (but has an exponential maximum on the brane) [3]. In infinite bulk 
models KK continuum begins from m = and this affects the long range behavior of gravity and 
other massless fields [10] [19] [20]. Massive fields if they have bulk modes can decay to the bulk with 
a life time which depends on the fundamental scale of gravity [19] [20]. For orbifoldized models the 
spectrum of KK modes is discrete. The long range effect of massive graviton modes is less important 
but the probability of decay of massive modes to the bulk is unchanged (see next section). Universal 
extra-dimension models in which all SM particles propagate to the compact dimensions are not ruled 
out for compactification radius of order TeV~^ or even lower [21]. 

At present brane world models have been constrained only based on the probability of direct obser- 
vation of processes involving the production of gravitons and its Kaluza-Klein modes [22] [23] [24] . A 
detail investigation of observable signal of the RS type models in Tevatron and LHC are performed 
in [25] and KK-mode production in the early Universe in [26]. The existent and near future accelera- 
tors can constrain the scale of gravity (and thus the size of the extra-dimensions) up to ~ SOTeV. The 
interaction of Ultra High Energy Cosmic Rays (UHECR) with protons in the terrestrial atmosphere 
has a CM energy close to lOOOTeF and is the most energetic interaction of elementary particles we 
can study today. It can be used to constrain the fundamental scale of gravity and the compactification 
scale up to much higher energies. 

The mass of KK modes detected by an observer on the brane is the result of smeared dimensions in the 
wave function. Classically however, it can be interpreted as a delay in the displacement of particles. 
For the observer on the brane if the delay slightly modifies the propagation of the particle in the 
detector, it is interpreted as a larger particle mass, otherwise it is seen as an arrival delay, specially 
with respect to the particles which propagate only on the brane. In the case of an air shower, the time 
coherence of the showers will be destroyed. 

In this work we calculate the minimum propagation time of particles ejected to extra-dimensions 
for a number of warped brane models and compare them with arrival time resolution of present Air 
Shower detectors. We restrict our study to = 5 models. This is enough for understanding the 
general characteristics of the propagation in extra-dimension from the point of view of an observer on 
a (3 -|- l)-brane and can be considered as a special configuration for models with higher dimensions. 
Our attention is mostly concentrated on the classical structure of the brane models because results are 
independent of the detail of quantum field contents and origin of the branes. However, before doing 
this we must assess the possibility and the probability that UHECRs' interaction in the atmosphere 
can produce bulk modes i.e. escaping particles. Given the rarity of UHECRs, only models in which 
the probability of production of bulk modes is very high can be constrained by this method. 

2 Production of Bulk Modes 

The whole idea of constraining brane models with UHECRs depends on the possibility and probability 
that remnants of UHECR interaction in the atmosphere can penetrate into the extra-dimensions. In 
this section we review the localization of fields on the brane with a special attention on models which 
provide bulk modes i.e. not all fields are confined to the visible brane at all energies. We estimate the 
probability of producing these modes at energy scale of interaction between UHECRs and nucleons in 
the terrestrial atmosphere. 

By definition in universal models [21] any field has bulk modes and propagates in all space-time 
dimensions. The interesting case for solving the hierarchy problem is when the size of the compactified 



dimensions are of the order of weak interaction or lower . This is ~ 3 orders of magnitude smaller 
than the CM energy of UHECRs' interaction. Therefore in the case of universal models, UHECRs can 
produce lowest KK-modes abundantly. In [23] it is argued that UHECRs can not probe the physics 
at very high energies simply because their interactions is dominantly electromagnetic. It is true that 
probability of the exchange of a heavy particle e.g a massive boson related to symmetries beyond 
Standard Model is very small with respect to a low transverse momentum EM cascade. However, 
in universal models as SM fields have bulk modes at very high energies all dimensions are "seen" 
as to be the same and the parameter space of EM cascades with non-zero momentum component in 
the extra-dimensions is much larger than cascades restricted to the three infinite dimensions. In the 
language of KK-modes, with a good accuracy the lowest modes can be considered as massless and 
they can be produced with the same probability as zero modes. In the following subsections we argue 
that for non-universal models one expects that at the CM energy of UHECRs interaction, most of 
localization mechanisms be no longer active and Standard Model particles can escape to the bulk. 



2.1 Scalars and Spin-2 Fields 

In warped models scalars and spin-2 fields can be confined geometrically. The zero mode of these 
fields has an exponentially decreasing wave function in the bulk [3] [27]. For models with infinite bulk 
the KK-spectrum begins from zero and therefore there is not a real confined zero mode. In addition, 
if the field has a non-zero 5-dim. mass, it has been shown [11] that 4-dim mass eigen modes on the 
brane are complex and decay to the bulk with a width: 

r/m4 oc (771,4/^)^ m\ = m\/2 (1) 

where 777.4 is the real part of the 4-dim. mass eigen value, 7775 is the 5-dim. mass of the field, and ^ is 
the warp scale in the static RS metric (Eq. (|56j) below). If the fundamental scale of Quantum Gravity 
is comparable to the weak interaction scale, at CM energy of UHECRs interaction it is expected that 
due to radiative corrections even massless particles like gravitons have an effective non-zero mass. For 
short distances relevant to the propagation of UHECRs in the terrestrial atmosphere, massive modes 
have a Yukawa type potential and their coupling is exponentially suppressed [11]. However, if the 777,4 
is much smaller than CM energy, the effect of exponential term in Yukawa potential is negligible. For 
most energetic UHECRs Ecm ~ IQ^^eV . This means that the coupling to modes as massive as ITeV 
is roughly the same as massless modes. The width in depends on the effective 5-dim. mass of the 
field and the warping scale. We discuss their implication on the decay of massive modes to the bulk 
and on the test of brane models in Section HI 

The above argument is also true in the case of 2-brane models where the spectrum is discrete. We 
show briefly that in static/quasi-static models the zero mode of scalar/spin-2 fields with 7775 > 
has an imaginary part i.e. it decays to the bulk (See also [28]). We determine the propagator of 
scalar /graviton using the Green function method discussed in detail in [27]. 

After changing variables y in metric (|56j) to: 

z = -e^y (2) 

^J' 

ds^ = -^(r]^^dx^'dx'' - dz^). (3) 

we apply the boundary conditions to both branes. Without loss of generality we assume one of them 
is at y = or 2; = ^ = i? and the other at y = L or z = ^e^^ = R' ^. The Green function (2-point 

^In some of universal models the extra-dimensions are not warped. Here we only study the propagation in warped 
spaces. Nevertheless, when the compactification scale is much larger than fi (see Sec. for definition), the warp 
factor is very close to one and the results of following sections are applicable. 

^We are only interested in the case where bulk is static. Therefore implicitly it is assumed that these points correspond 
to the fixed points of the radion field. 



propagator) A{x, z,x' , z') is the solution of 4-dim. mass eigenstate equation: 



{z'^d1 + zd,+p'^z^-d^)Ap{z,z') = t^5{z-z') (4) 
A{x,z,x',z') = l^e'P(^-^'^A,{z,z') (5) 

^ (^) Ap(z,/) (6) 



d = ^A + R^ml (7) 

The boundary and matching conditions for right and left propagators: 

A< = Ap{z, z') z < z' , A> = Ap{z, z') z > z' (8) 

are as followings (boundary conditions are deduced from Q and matching condition from (jU): 

d,iz^A>)U=ji, = (9) 
d,{z^A<)\,=R = (10) 
A<|^=^/ = A>\,=,> (11) 

5,(A<-A>)U=,, = ^ (12) 

Solutions of @ are linear combination of Bessel Functions: 

A{z', R, R')Jd{pz) + B{z\ R, R')Nd{pz) (13) 

Applying the conditions (|^- (|12l) to this solution leads to an equation which determines the KK mass 
spectrum: 

pR' ,UpR') + {I - u)J,+i{pR') ^ pRUpR) + {I - y)J,+i{pR) 
pR'N^{pR') + {l-u)N^+i{pR') pRNuipR) + {I - v)N^+i{pR) ^ ' 

with V = d — 1. Finding the exact solution of (|14() is not trivial. To consider only a simple case we 
assume that 4-dim. mass of the scalar field \p\ <^ ^ = i.e. pR <C 1. Regarding the Standard Model, 
this can be applied to a confined Higgs when the scale of compactification is much higher than Higgs 
mass or to a light axion like scalar or to graviton with a small mass due to radiative corrections. We 
keep only lowest powers of pR in the expansion of J^,. For solving hierarchy problem R' ^ R. Using 
the asymptotic expansion of Bessel functions, (fTI|) reduces to: 



. (15) 

zy2-''sin(i/7r)r(-i/) ^ ^ 

where t] = pR and rj' = pR' . For ^ 1 the solutions of ()15() lead to the following mass spectrum: 

ipo = m4 — iT (16) 

m4 « fiy^2u{u-l) « ^ , rp^msVs^ (17) 

v2 

w ^ (n7r+ — + — ) For large ?i. (18) 

Continuity properties of guarantees that H15|) is also valid when ^ or equivalently z/ — > 1. 
In this case = or = 0. For pR <C 1 the mass difference between KK-modes p„ is Ap oc 
1/R' <ti 1/R = jJL. Due to special properties of Bessel Function with integer index, the zero mode of 
massless fields is protected from decay even when the spectrum of KK-modes for massive particles 



begins roughly from zero. Tunneling probability depends on 5-dim. mass of the field and on warping 
scale ^. If /i is large, the probability of zero-mode decay to the bulk can be small. However, except for 
very light particles a /i as large as ITeV provides an enough large width (F > W~^^eV) for decay to 
the bulk during propagation in the terrestrial atmosphere if is in the mass range of SM particles. 

To see the effect of mass on the coupling we can investigate the mass dependence of the propagator 
on the branes. Using (jH)) and ©-dJ, one can determine the integration coefficients A{z' , R, R') and 
B(z' , R, R') in and right and left propagators: 



A<(z,z') 



A>(z,z') 

Ao 
Ai 
A2 
A3 



7rR:^(^AoJ^+i{pz') - AiN^+iipz')^ (^A2Ju+iipz) - AsN^+iipz)^ 



z'2(AoA3 - A1A2) 



7rR^(^A2Ju+i{pz') - A^N^+iipz')^ (^AoJ^+i{pz) - AiN^+i{pz))^ 



^'^(AoAa 
pR'N^{pR') + (1 - v)N^+i{pR') 
pR'MpR') + (1 - u)J^+i{pR') 
pRN^ipR) + (1 - u)N^+i{pR) 
pRMpR) + {1 - u)J^+i{pR) 



A1A2 



Restricting these equations to the branes gives the 2-point propagators: 



A{x,R,x',R) 
Aix,R',x',R') 



1 



(27r)4 
R'^ 
R^{2^Y 



,_^,.p-\AQjy+i{pR) - AiN^+i{pR)) 



dp e 



(A0A3 - A1A2) 
Aip{x-x') P~\A2J,+i{pR') - A-iK+^{pR')) 



(19) 



(20) 

(21) 
(22) 
(23) 
(24) 



(25) 
(26) 



(A0A3 - A1A2) 

The term R!'^ jR? in ()26() reflects the difference between metric on the branes. Consistently, the roots 
of dominator in H25() and H26|) . are the same as H14|) and correspond to KK-modes. Near each mass 
mode the propagators can be written as a Yukawa propagator with complex mass and a coupling equal 
to the residue of the integrand. Applying this procedure to (|26|) one finds that on the brane at i?' for 
v>\: 



9l 
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(27) 



In our approximation \po/fJ.\ oc m^/fi <^ 1, \p'^/ fi\ > |po//^| and it can be even larger than 1. Therefore 
coupling to KK-modes is larger than to zero-mode. The probability of production of zero-mode with 
respect to n''^ KK-mode is: 

2(j/-l) / 2 

Y ^ PO Pji (28) 

Pn 

For KK-modes with \p'^\ ~ /_f the branching ratio V < 1 and therefore the probability of production 
of these modes is larger than zero-mode. 

In conclusion, radiative corrections that can induce an effective mass for scalar/spin-2 fields weakens 
their confinement on the brane. In warped 2-brane models even when 7715 = the mass difference 
between the zero mode and massive KK-modes on the brane with smallest warp factor is small and 
they can be abundantly produced in high energy interactions. 



2.2 Fermions and Gauge Fields 



Fermions can not be localized to the TeV scale brane (i.e. brane at R') gravitationally but a chiral 
symmetry breaking can confine them [14]. Their escape to the bulk when their 1715 > has been 



studied in detail for one brane models in [11] and [29]. The detail of the formalism is very similar 
to the case of a scalar field except that the mass eigenstate equation includes a chiral mass term due 
to interaction of fermion modes with a bulk scalar responsible for the symmetry breaking. The width 
of the zero mode depends on the coupling between fermions and the inferred scalar field and without 
detail knowledge of underlying particle physics it is difficult to assess the probability of decay to the 
bulk. For 2-brane models the situation should not be very different and general conclusions of Sec l2.1l 
must be applicable. 

At present the particle physics models don't fix the scale of the symmetry breaking. For not creating 
a new hierarchy however it can not be much larger than compactification scale [18] or fundamental 
scale of gravity M5. Therefore, at energy scale of UHECRs interaction in the atmosphere, not only 
it is possible to produce KK-modes, it is very probable that at such energies the restoration of chiral 
symmetry completely removes the confinement of fermions and open the extra-dimension even to 
fermionic zero modes (presumably SM matter). 

As for gauge bosons, the most successful scenario for their confinement on the branes is based on 
adding an induced kinetic term to the action of bulk gauge fields on the brane. It appears due to 
the interaction of these fields with confined charged scalar or fermions on the branes [17]. Other 
suggestions are mostly equivalent to this scenario [30]. Once the charged fields become able to escape 
to the bulk, they drag their interaction vertex with gauge fields to the bulk and release them from 
confinement. It has been shown [31] that the coupling of gauge boson KK-modes to fermions on the 
brane is stronger than coupling of their corresponding zero-mode (similar to the self coupling of scalars 

The general conclusion of this section is that regarding: 

- Very high CM energy of interaction of most energetic Cosmic Rays in the atmosphere which is 

much higher than natural confinement scale of Standard Model particles on the brane and natural 
fundamental scale of gravity to solve the hierarchy problem; 

- The fact that confinement of SM fields is not intrinsic but the result of either a broken symmetry 

(for fermions) or interactions (for bosons); 

- That the particle physics in 5- dim can not be completely massless and at least part of the particle 

spectrum must acquire mass as it is the case in observable 4-dim Universe. In fact as radion 
is in fact the scalar component of 5- dim. metric perturbations [9], it couples to all bulk fields 
and radiatively induces a small mass term. Consequently, zero-modes (presumably SM particles) 
are not stable and in a finite time decay to the bulk unless another phenomenon like symmetry 
breaking prevent it; 

the phase space of production of bulk modes at high energy tail of UHECRs spectra seems to be higher 
than confined modes. 

Until now more than 100 coherent air showers have been observed with Ecm ^ SOOTeV (assuming 
interaction with nucleons in the atmosphere), 17 between them have energies more than 450Tey and 
one has a CM energy close to lOOOTey [32]. Assuming that UHECRs interaction with Ecm ^ 300Tey 
produces bulk modes abundantly, the mere observation of coherent showers up to energies close to 
10'^^ eV constraints the parameter space of brane models. 

As the assessment of cross-sections and other details are model dependent, in the rest of this work 
we simplify the problem of constraining brane models and consider only the classical propagation of 
the bulk modes. This method has been already used by other authors to study some of cosmological 
consequences of brane models [33] [10]. 



3 Propagation 



The geodesic path of particles in the bulk has been already studied in a number of previous works [33] [10] [11] 
[34]. However, most of them are concerned with the possible acausality of paths for observer on a 
brane and their purpose is to see if it can solve the cosmological horizon problem in the early universe. 
Here we are concerned with the present evolution of the Universe and simplifying assumptions will be 
based on its present very slowly changing state. 

The metric of the 5-dim brane models can be written as the following: 

ds^ = n^{t, y)dt^ - a^{t, y)5ijdx'dx^ - b^{t, y)dy'^. (29) 

For a static bulk b'^[t,y) is constant and we can normalize coordinates such that 6=1. The geodesic 
of a particle is defined by: 



duP n n n 2n' 



+ + uV + ^nV5i^- = (30) 



aa 



dr ' n ' n ' 
du^ 2a' ,• A 
dr a a 



du" 2a' i ^ 2d i Q 

— u'u^ H u'u^ = (31) 



dii* 

+ W-uV - aa'uV<5i,- = (32) 

dr 

dt i dx' 4 dy 

= u = —, u =-r- (33) 

dr dr dr 

T is the proper time parameter along the particle world line. It is easy to see that H31|) is integrable 
and: 

u' = % (34) 

where ^* is an integration constant. As we are only interested in the minimum delay in the arrival 
of particles due to the propagation in the extra-dimensions, we put ^* = 0. Later we try to estimate 
qualitatively the effect of a non-zero . Even after this simplification the system of equations pOj) - 
(|32jl is highly non-linear and coupled. In the following we calculate an analytical solution for the case 
njn ~ 0. This approximation is justified when we are interested in the propagation of particles in an 
extremely short period of time with respect to the expansion rate of the bulk or the brane. In fact 
from the solution of the Einstein equations [35], n{t,y) can be normalized such that: 

nit,y) = —j-^ (35) 
ao{t) 

where ao{t) = a{t,y = 0) assuming that one of the branes is at y = 0. At present, both a(t, y) and 
do(t) are very slowly varying quantities. Therefore n{t,y) ~ n{t(),y) where to is the present time, and 
h{t,y) ~ 0. 

Under this approximation: 



and: 



^ « n'u^ (36) 
dr 



du^ 2u° dn , , 

— + — — = 37 
dr n dr 



du'^ n{v?)'^ dn 
dr dr 



V\ Lb Lit l> - > 

+ ^- = (38) 



The solution of (|H7|l and is straightforward: 

u° = en-'^, = ±{i] - ^] (39) 



The parameters 9 and rj are integration constants and must be determined from the initial conditions. 
The lb sign defines the direction of the propagation. In the rest of this letter we neglect the direction 
and consider only the absolute value of u'^ . 

For a particle leaving the visible brane placed at y = yi, at time t = to: 

q2 Massive particles, 

9 = u^ito,yby{to,yb), ^ = ^77 7 - (^^^(^0, y?,))^ = <^ (40) 

IT' v-o^yb) LO Massless particles. 

After eliminating the proper time from and u^, one obtains the equation of motion in the bulk (For 
simplicity we assume that Dy/dr k, dy/dr and Dt/dr « dt/dr): 





= n\t,y)alit) = 


Ait) 


= ao\t)-Cit), 


Bit) 


= -Pfeo«o^(i), 


Cit) 


2d2(t) 

2 ' 












= ]I—\pb\ 



dy _ n\t,y) / 
dt 9 ^n^it,y) 

•1 Massive particles, 

(42) 

-0 Massless particles. 

Our approximations are valid only when dy/dt is real. This put limits on the testable part of the 
parameter space of the models (see below). 

Einstein equations give the solution for a(t,y) and n(t, y) [35]. For a flat visible brane: 

a^it,y) = ^(t)cosh(/iy) + e(t)sinh(/xy) +C(t), (43) 

Ait) cosh(/iy) + Bit) sinhii^iy) + C(t) 
4a2(t,y) 

(45) 
(46) 

(47) 
(48) 

For any density p, p' = p/kus, Ars = •ip/K?- The densities and pB are effective total energy 
density of the brane at y = and the bulk respectively. We consider only AdS bulk models with 
Pb < 0. The constant k'^ = Stt/M^ is the gravitational coupling in the 5-dim. space-time. The model 
dependent details like how p'^^ and pB are related to the field contents in the bulk and on the brane 
and how they evolve are irrelevant for us as long as we assume a quasi-static model. The solution 
(|43|) is valid both for one brane and multi-brane models. The only difference between them is in the 
application of Israel junction conditions [15] [16]. 

Equation (|42|) is non-linear and its integration non-trivial. We use again the quasi-static properties 
of the present Universe and its low energy density to simplify the integration. Matter density on the 
branes at late time is much smaller than the brane tension or induced tension by scalar fields [38]. 
Therefore, it is not unreasonable to neglect time dependence of densities in H45|) - (|47|) and to consider 
only cosmological constant type energy-momentum densities. This simplification is even more justified 
in our case where we have to deal only with very short duration of the propagation in the extra- 
dimensions. This approximation and (|45j) - (|47j) lead to: 

C = -i?^. (49) 

A = 2ao(t)do(t) - Cit), (50) 
B = -2p[^aoit)doit). (51) 



After changing variable y to z = e^^ and using H35|) : 



a^{t,z) 



V = - 



2z 
(1 - P^o 



V{t,z), 



C{t). 



z^ + 



2C{t) 



z + + p'b. 



C{t). 



dz 
IE 



(52) 
(53) 

(54) 



If an ejected particle to the bulk comes back to the brane, must go to zero at some point in the 
bulk before the particle arrives to the bulk horizon (if it is present). The roots of 1)54(1 correspond to 
these turning points and determine the propagation time in the bulk. In the next simplifying step we 
use again the fact that the typical propagation time we are interested in is very much shorter than 
the age of the Universe and therefore B, C and A, B, C during propagation are roughly constant, the 
right hand side of ()54() depends only on z and is easily integrable: 



At. 



propag 



2{t 



stop 



to) 



Zstop 



2dz 



(55) 



In ()55|) . to is the initial time of propagation in the extra-dimension and tgtop is the time when the 
particle's velocity changes its direction, i.e. when dz/dt = 0. The integral in ()55() is related to the 
elliptical integrals of the first type J-{iv, v) where lo and v are analytical functions of the denominator 
roots in and zq [39]. Note that Zgtop corresponds to the closest root to zq. 



4 Test of Brane Models 

In this section we apply the formalism discussed in the previous section to most popular brane models 
and determine the propagation time of high energy particles in the fifth dimension. Note that the 
calculation of propagation time in these models under our approximations is valid only for durations 
very smaller than the age of the Universe and if in the following figures in part of the parameter space 
the propagation time can be larger, this part of the figure should not be considered. 



4.1 2-Brane Models 

It has been shown in [15] [16] that by imposing constraints on the visible brane to obtain the observed 
value of cosmological constant A and Newton coupling constant G and to solve the hierarchy problem, 
all parameters of this class of models i.e. poj PL and Ajj5, can be determined as a function of 
where L is the distance between two branes. It is not however possible to find an exact analytical 
form for the solutions. Moreover, the analytical solution in [16] has been obtained for a special setup 
which decouples hidden and visible branes. Here we free some of the constraints, first to be able to 
find analytical solutions, and second to extend this study to a larger number of models. 



4.1.1 Geodesies in Static RS Models 

In the original RS model with static metric: 

ds^ = e-^'yi^^^dx^'dx" - dy\ (56) 

the main constraint on the model is the cancellation of the cosmological constant on the visible brane 
which leads to the equal and opposite sign tensions po = —pi = Ars- The solution of the hierarchy 
problem limits the range of the parameter fiL. It has been shown [10] that in the fine-tuned RS 
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Figure 1: Propagation time for relativistic particles with u^{to)/N = 10'^ (full line) and u^{tQ)/N = 
1.2 (dash line) in RS model. Red, magenta and light green curves correspond to M5 = lO^^e^, 
M5 = lO^^eV and M5 = lO^^eV (or fi ~ IQ-^^eF, fi ~ IQ-^eF and fi ~ W'^eV for fine-tuned 
model) respectively. The dark green line shows the time coherence precision of present Air Shower 
detectors. 



model photons leave the brane and never return. Assuming the visibility of the all dimensions of the 
space-time at very high energies, in this model we could never observe the ultra high energy particles 
and therefore it is automatically ruled out. Nonetheless, to test the formalism of the previous section 
we apply it to this model. 

For a very small cosmological constant (as it is assumed in the process of fine-tuning) C(to) ~ and 
P(to,2;) ~ 1. For massive particles, the denominator in (|55j) has only one root: zi = 1/6 and: 



= ^f-(^ ^ = e-M(^o). (57) 

Fig^ shows Atpropag as a function of fiL and A/5 for massive relativistic particles. With present air 
shower detectors time resolution of order W^sec, only when M5 > lO^^ey, the model is compatible 
with the observed time coherence of the UHE showers. For fine-tuned RS model fi ~ G/k^ [3] i.e. 
II = GM^ ~ W~^eV for M5 ~ 10^^ eV [3]. Due to smallness of ji and consequently lightness of 
KK-modes for SM particles even this model with large M5 has already been ruled out [31] unless a 
conserved quantum number prevents the production of KK-modes [21]. 

For massless particles: 

z{t) -zo = fi{t - tof (58) 

In (|58|) . z{t) is monotonically increasing and there is no stopping point. With our approximations 
there is no horizon in the bulk because a{t, y) is roughly constant. Therefore (|58() means that massless 
particles simply continue their path to the hidden brane and their faith depends on what happen to 
them there. At very high CM energy of UHECR interactions if charged particles can escape to the 
bulk photons are also dragged to the bulk and never come back. 



4.1.2 General Solution 

Numerical solution of constrained 2-brane models in [15] [16] shows that for /iL > 5 the tension on 
both branes is positive and very close to A/jg. We can use constraints on the Cosmological Constant 
and hierarchy to find pg and = p'^. We redefine them as pg = 1 + Apg and p'^ = 1 + Ap'^. To solve 



hierarchy problem (See equations 29-31 in [16]): 

Mlr^N^='A= PAo(l-cosh(/iL))+sinh(/iL) ^ ^ 
M^l^ nl p^^(l-cosh(^L))+sinh(^L) 

This leads to: 

1 _ e-Mi + Ap'^(l - cosh(/.L)) ^ _ 



iV 



" 1 - cosh(^L) 1 - cosh(^L) 

For a very small A^^ and Ap^ < 1, the first term in ()6U() is ©(A^^) and: 

"l-cosh(/iL) "l-cosh^^L) ^^^^ 



Using a\/a?j^ = where H is the Hubble Constant on the visible brane [16]: 



2iV2 sinh(/iL) 



(1 - cosh(^L))^5- + e~^^ - 1 ± 



/ 2i/2 2iJ2 

((1 - cosh(//L))^5- + e-^'^ - if + iV^ sinh(/iL)(^5- + 26"^^^) 



(62) 



In (|62|) the solution with plus sign gives Ap^ ~ —2 which deviates from our first assumption |A/j^| < 1 
and leads to a negative tension on the visible brane like static RS model. The solution with negative 
sign is: 

^p'l ^ 2e-f'^ (63) 
and both branes have positive tension close to Ajis- 

When the matter densities on the branes and in the bulk are negligible [16], Oq/^o ~ ^l/'^'l ~ 
and: 

,64) 

It is easy to see that V and a?{t,y) have the same roots. Models with a horizon i.e. a point such 
that a?{t,yfi) = are pathological (because no particles/brane behind it is observable). The condition 
to have no real root i.e no horizon in the bulk is: 

-2 + — <Ap^ + C'^-— (65) 

For massive particles, the denominator of the integrand in (|55|) can have two roots: 

_ C-Q^± ^{C - ff^y + (2 + Ap(, + CO(Ap^ + CO 

Ap^ + C ^''''^ 

The model is consistent only if I?(z — P/0^) > in the range of integration. Therefore: 

z+ ^ zo = e''-^ ^ 2-- (67) 

The matter on the brane is confined only if > 1 and Zj^ zq when — > 0. To first order in e~^^^ 
and A^^ this leads to the following relation between parameters of the model: 

-C' = ^p', + -{N^ + ^p',) (68) 

This condition is not an addition to the model described in [16]. It is in fact the result of solutions 
H61|l and (|63p under the approximations considered here. It is not evident whether such a constraint 
appear in the full theory. 
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Figure 2: Left: Propagation time for relativistic particles in 2-brane model. Description of the curves 
is the same as Fig^ Right: Parameter n{eV) as a function of fiL. It is roughly independent of M5 
(See the text). 



For fiL > 5 the right hand side of (|H5|) is positive. Therefore —C oc H"^ can not be zero. This 
relation between a small but non-zero value of the Hubble Constant or equivalently Cosmological 
Constant on the visible brane and the smallness of N'^ and fi which is related to the strength of the 
induced gravitational coupling on the brane, confirms the same observations in [16] for an analytical 
solution of 2-brane models with some approximations and in [37] for the exact solution of some special 
models. 

Finally the propagation time in the bulk is given by: 

4 

^tpropag = -T{a,Q). (69) 

^,mp', + C'\)^ 

a = 2arctan (70) 



Q 




C'{C' - 02) + 4(2 + Apf, + CO(Ap(, + C) 



2 / ^/ x2 



(71) 



Fig m shows the propagation time for models which satisfy simultaneously (|61() . (|63() and H68|) . In 
equation (|68j) up to first order, C' depends only on fiL and thus in (|64j) the value of is independent 
of M5. Only models with large fiL ^150 are not ruled out. This is due to (|65j) and smallness of the 
observed Hubble Constant H. The same conditions make fi > leV which is much higher than the p 
for the fine-tuned RS model. With // < rriradion [6] [9] these class of brane models are also consistent 
with constraint on the fifth-force measurements [41]. From (|17() the corresponding life-time for SM 
fields with 7715 > is shorter than W~^^sec, much shorter than propagation time in the atmosphere 
and also much shorter than propagation time in the bulk. This justifies the classical treatment of 
propagation. 



We have also applied the formalism described here to the fine tuned model of [16]. In this model 
the equation of state on the hidden brane is fine-tuned to neutralize its effect on the visible brane. 
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Figure 3: Left: Propagation time for relativistic particles in the fine-tuned 2-brane model of [16]. 
Right: Parameter /i as a function of M^. 



This results to a unique definition of G and the only free parameter in the model is M5. Although 
this model has been obtained just by phenomenological arguments, a field theory model suggested 
by Arkani-Hamed et al. [36] to solve the Cosmological Constant problem has the same form for 
G if iiL ~ k?(f){Q) where (/)(0) is the vev of radion on the visible brane (The Arkani-Hamed et al. 
model has only one brane but it includes a horizon in the bulk which limits the accessible size of the 
extra-dimension and makes it similar to a two brane model). 

FigElshows the propagation time and /u as a function of M5. Unfortunately despite physical interest of 
this model it is only compatible with very high M5 > 10^^ el^ unless the decay to the bulk is prevented 
up to such high energies. The corresponding however is consistent with present constraint on the 
Fifth force [41] and the life-time of zero-modes of massive 5-dim. fields for such high M5 is enough 
short to permit particles to decay to the bulk during their propagation in the atmosphere. The value 
for M5 is some 6 orders of magnitude larger than lOTeF, presumably the Electroweak interaction 
scale and it would be a matter of speculation to consider this model as having no hierarchy problem. 
Another problem in testing these models with UHECRs is that as the natural scale of gravity M5 is 
very high, it is possible that the symmetry breaking scale which is necessary for the localization of 
fermions (and indirectly gauge bosons) is also much higher than CM energy of UHECRs interaction. 
In this case only very weakly interacting particles like gravitons can decay to the bulk. As the total 
production cross-section for them can be tiny, the number of observed UHECRs event can be not 
enough to constrain such models. We have also tested the general 2-brane models without taking into 
account (|68)) . Roughly speaking, it is equivalent to having a comparable matter density and tension on 
the hidden brane. The value of jj. becomes a free parameter. The result is shown in Fig0]for 3 different 
values of /i. Models with fi > IQ^eV and /iL < 70 are compatible with the present observation of 
UHECRs. The lower limit for ^ from this test is higher than the constraint obtained from Fifth force 
experiments [41]. 

4.2 One-Brane Models 

The solution of Einstein equations for symmetric one-brane models is the same as two-brane ones [35] . 
Due to existence of only one boundary however the bulk and brane tensions are not related. In 
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Figure 4: Propagation time for relativistic particles in 2-brane models with as a free parameter. 
Top left: fj, = 10~'^eV; Top right: /x = lO^eF; Bottom: fi = 10^ eV. Description of the curves is the 
same as FigQ] 



addition, the cosmological evolution on the brane: 



-^PB + — {pb + Pm[t)) + 

3d 



C{t) 



(73) 



includes an arbitrary function C{t) which is related to the bulk tension and matter [40]. Here we test 
two popular models studied in [35] and [40]. 

In the first model [35] C{t) = and the brane tension is fine-tuned to cancel the effect of quadratic 
term at late times: 



^2 



and: 



SvrG 



-4 

K Pb 



(74) 



(75) 



The only free parameter in the model is M5. In the second model [40] C{t) (or equivalently T5 
component of the energy-momentum tensor) is adjusted such that the conventional evolution equation 
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Figure 5: Left: Propagation time for one-brane models of [35]. Right: Parameter /x as a function of 
fiL. 



be obtained. At late times when the brane tension is much larger than time dependent matter terms 
these two models are roughly the same. 

Equations ((71]) and ((75|) determine pf, and ^. It is easy to see that pf, = A^^. The definition of C 
and roots are the same with Apo = ^Pb = and zq = 1. FigEl shows the propagation time for these 
models. 



4.2.1 Effect of 

When 0* / equations H37|) and can be written as: 



+ 



dr n dr 
, du^ a'O'e^S, 



+ 



dn 



+ {vl'Yn— = 
dr 

At least formally Ea. H76|l can be solved analytically: 

_ e ^ e'e^,jg{t,y) 



u 



+ 



dr[-^ 



(76) 
(77) 

(78) 
(79) 



We don't need to solve 1)77(1 directly. Knowing vP and we can use the definition of velocity vector 
to determine u^: 



n 



2(«Y-aVu^5,,-(u^)^ 



JlJ \^Ub J S (^^) 

The definition of e is the same as in (|42() . After elimination of dr we obtain the formal description of 
the equation of motion in the bulk: 



dt 



tj \ ' f^Z 



a- 



, QB^ei6,j 

— 7 ~r ^ 



(81) 



We can use (jHl]) to determine dr (As before for simplicity we assume that Dx'^/dr « dx^/dr): 



and: 

Q'e^bifi{t(T),y(T)) = -J ^6,j6^dx' (83) 

In d/a is independent of x*. The rest of right hand side of (|S5|) i.e. / 6ij6^dx^ is the projection of 
the particles world line on the brane. The value of Hubble constant d{t, y)/a(t, y) ~ (from (|45I47|) 
and H49I51|) ') and thus 6^9^5ijQ is very small when the projection distance traversed by the particle 
is small with respect to the Hubble radius. Therefore we presume that conclusions of the previous 
section will not be extremely modified when full propagation is considered. 



5 Conclusion 

The calculation in this work is mainly based on two assumptions: 

- At interaction energy scale of most energetic Cosmic Rays the physics is high dimensional either 

because all symmetry based confinements are no longer at work or because there is a large 
cross-section and/or phase space volume which permits the production of massive KK-modes. 

- In the time scale of the propagation of a particle in the extra-dimension, the bulk and the branes 

are quasi-static. 

If these assumptions are valid the time coherence of Ultra High Energy Air Showers rules out a large 
part of the parameter space for a number of brane models unless some micro-physics phenomena 
confine particles to the brane at energies much higher than Electroweak scale. This makes a new 
hierarchy inconsistent with the spirit of brane models. 

For most 2-brane models the acceptable range of ;U is /i > \eV except for original RS model which 
needs n > 10~'^eV. This lower limit is much lower than what can be obtained from non-observation of 
KK-mode production in accelerators [31]. However, at present accelerator energies it is always arguable 
that presence of some quantum conservations prevent the production of KK-modes. Presence of such 
conservations at the CM energy of UHECR interaction seems much less natural and constraints on 
M5 more are robust. 

The upper limit of L is also a universal value for models with different range of fiL: L < 10~^eV~^ 
10~^ cm. It is much smaller than the upper limits obtained from gravity experiments [42] [41]. Again 
for static RS model the upper limit is L < lO^cm, much larger than one obtained for other models. 
One brane models with interesting range of M5 are ruled out. 

This study has an additional interesting conclusion: The close relation between a very small but non- 
zero Cosmological Constant and the smallness of the Newton coupling constant (i.e. the hierarchy 
problem). In fact without the fine-tuning of these apparently independent physical quantities, the 
brane models are not consistent. 
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